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C/^ ! Abstract 

Asymptotic stabihty of smaU sohtons in one dimension is proved in the framework of 
a discrete nonhnear Schrodinger equation with septic and higher power-law nonhnearities 
and an external potential supporting a simple isolated eigenvalue. The analysis relies on the 
dispersive decay estimates from Pelinovsky & Stefanov (2008) and the arguments of Mizu- 
machi (2008) for a continuous nonlinear Schrodinger equation in one dimension. Numerical 
^ , simulations suggest that the actual decay rate of perturbations near the asymptotically 

■ stable solitons is higher than the one used in the analysis. 



1 Introduction 



Asymptotic stability of solitary waves in the context of continuous nonlinear Schrodinger equa- 
tions in one, two, and three spatial dimensions was considered in a number of recent works (see 
Cuccagna [J| for a review of literature). Little is known, however, about asymptotic stability 
of solitary waves in the context of discrete nonlinear Schrodinger (DNLS) equations. 
^ ' Orbital stability of a global energy minimizer under a fixed mass constraint was proved by 

^ . Weinstein [24] for the DNLS equation with power nonlinearity 

iUn + AdUn + \Un\'^^Un = 0, n ^ if" , 

where is a discrete Laplacian in d dimensions and p > 0. For p < ^ (subcritical case), it 
is proved that the ground state of an arbitrary energy exists, whereas for p > | (critical and 
supercritical cases), there is an energy threshold, below which the ground state does not exist. 

Ground states of the DNLS equation with power-law nonlinearity correspond to single- 
humped solitons, which are excited in numerical and physical experiments by a single-site initial 
data with sufficiently large amplitude [IT]. Such experiments have been physically realized in 
optical settings with both focusing [7j and defocusing [13] nonhnearities. We would like to 
consider long-time dynamics of the ground states and prove their asymptotic stability under 
some assumptions on the spectrum of the linearized DNLS equation. From the beginning, 
we would like to work in the space of one spatial dimension [d = 1) and to add an external 
potential V to the DNLS equation. These specifications are motivated by physical applications 
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(see, e.g., the recent work of [T6] and references therein for a relevant discussion). We hence 
write the main model in the form 

iUn = {-A + Vn)Un + 'y\Un\'^^Un, 71 G Z, (1) 

where Aun ■= Un+i — 2un + Un-i and 7 = 1 (7 = —1) for defocusing (focusing) nonlinear- 
ity. Besides physical applications, the role of potential V in our work can be explained by 
looking at the differences between the recent works of Mizumachi [K] and Cuccagna [S] for 
a continuous nonlinear Schrodinger equation in one dimension. Using an external potential, 
Mizumachi proved asymptotic stability of small solitons bifurcating from the ground state of 
the Schrodinger operator Hq = —d"^ + V under some assumptions on the spectrum of Hq. He 
needed only spectral theory of the self-adjoint operator Hq in since spectral projections and 
small nonlinear terms were controlled in the corresponding norm. Pioneering works along the 
same lines are attributed to Soffer-Weinstein [20l EH |22] , Pillet & Wayne [19] , and Yao & Tsai 
[25} [26l [27] . Compared to this approach, Cuccagna proved asymptotic stability of nonlinear 
space-symmetric ground states in energy space of the continuous nonlinear Schrodinger equa- 
tion with V = 0. He had to invoke the spectral theory of non-self-adjoint operators arising in 
the linearization of the nonlinear Schrodinger equation at the ground state, following earlier 
works of Buslaev &: Perelman [H [2] , Buslaev & Sulem [3] , and Gang & Sigal O |9] . 

Since our work is novel in the context of the DNLS equation, we would like to simplify 
the spectral formalism and to focus on nonlinear analysis of asymptotic stability. This is the 
main reason why we work with small solitons bifurcating from the ground state of the discrete 
Schrodinger operator H = —A+V. We will make use of the dispersive decay estimates obtained 
recently for operator H by Stefanov &: Kevrekidis [23] (for V = 0), Komech, Kopylova & Kunze 
|12j (for compact V), and Pelinovsky & Stefanov [18J (for decaying V). With more efforts 
and more elaborate analysis, our results can be generalized to large solitons with or without 
potential V under some restrictions on spectrum of the non-self-adjoint operator associated 
with linearization at the nonlinear ground state. 

From a technical point of view, many previous works on asymptotic stability of solitary 
waves in continuous nonlinear Schrodinger equations address critical and supercritical cases, 
which in (i = 1 corresponds to p > 2. Because the dispersive decay in — norm is slower for 
the DNLS equation, the critical power appears at p = 3 and the proof of asymptotic stability 
of discrete solitons can be developed for p > 3. The most interesting case of the cubic DNLS 
equation for p = 1 is excluded from our consideration. To prove asymptotic stability of discrete 
solitons for p > 3, we extend the pointwise dispersive decay estimates from [18] to Strichartz 
estimates, which allow us for a better control of the dispersive parts of the solution. The 
nonlinear analysis follows the steps in the proof of asymptotic stability of continuous solitons 
by Mizumachi [15] . 

In addition to analytical results, we also approximate time evolution of small solitons nu- 
merically in the DNLS equation ([T]) with p = 1,2,3. Not only we confirm the asymptotic 
stability of discrete solitons in all the cases but also we find that the actual decay rate of 
perturbations near the small soliton is faster than the one used in our analytical arguments. 

The article is organized as follows. The main result for p > 3 is formulated in Section 
2. Linear estimates are derived in Section 3. The proof of the main theorem is developed in 
Section 4. Numerical illustrations for p = 1,2,3 are discussed in Section 5. Appendix A gives 
proofs of technical formulas used in Section 3. 
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2 Preliminaries and the main result 

In what follows, we use bold-faced notations for vectors in discrete spaces ll and 1^ on Z defined 
by their norms 



Components of u are denoted by regular font, e.g. Un for n G Z. 

We shall make the following assumptions on the external potential V defined on the lattice 
Z and on the spectrum of the self-adjoint operator H = — A + V in P. 

(VI) V G ll^ for a fixed > |. 

(V2) V is generic in the sense that no solution ■j/>g of equation HipQ = exists in for 



(V3) V supports exactly one negative eigenvalue ujq < of H with an eigenvector tpQ £ P and 
no eigenvalues above 4. 

The first two assumptions (VI) and (V2) are needed for the dispersive decay estimates 
developed in |:.l8j. The last assumption (V3) is needed for existence of a family cf){uj) of real- 
valued decaying solutions of the stationary DNLS equation 



near to = ujq < 0. This is a standard local bifurcation of decaying solutions in a system of 
infinitely many algebraic equations (see [TJj for details). 

Lemma 1 (Local bifurcation of stationary solutions). Assume that V G and that H has an 
eigenvalue loq with a normalized eigenvector tj^Q G P such that \\tpQ\\i2 = 1. Let e := uj — loq, 7 = 
+1, and eo > be sufficiently small. For any e G (0, eg), there exists an e-independent constant 
C > such that the stationary DNLS equation (0) admits a solution 4>iuj) G C^([u;o!'^o + eo])^^) 
satisfying 




2 < (J < 2- 




(2) 




Moreover, the solution 4>{^) decays exponentially to zero as \n\ 00. 
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Remark 1. Because of the exponential decay of (f){Lv) as \n\ — > oo, the solution (f){uj) exists in 
l^ for all a > 0. In addition, since \\(t)\\ii < Co-H^Hia, for any a > ^, the solution ct){io) also 
exists in l^ . 

Remark 2. The case 7 = — 1 with the local bifurcation to the domain uj < ujq is absolutely 
analogous. For simplification, we shall develop analysis for 7 = +1 only. 

To work with solutions of the DNLS equation ([T|) for all t E M_|_ starting with some initial 
data at t = 0, we need global well-posedness of the Cauchy problem for ([T|). Because H is a 
bounded operator from P to P, global well-posedness for ([1]) follows from simple arguments 
based on the flux conservation equation 

^^I'UnP = Un{Un+l + Un-l) - + Un-l) (3) 

and the contraction mapping arguments (see [T7] for details). 

Lemma 2 (Global well-posedness). Fix a > 0. For any uq G l'^, there exists a unique solution 
u{t) G C"'^(M+,/^) such that u(0) = uq and u{t) depends continuously on uq. 

Remark 3. Global well-posedness holds also on M_ (and thus onM) since the DNLS equation 
(OP is a reversible dynamical system. We shall work in the positive time intervals only. 

Equipped with the results above, we decompose a solution to the DNLS equation ([T]) into 
a family of stationary solutions with time varying parameters and a radiation part using the 
substitution 

u(t) = e-^W(<^(..(t))+z(t)), (4) 

where (a;, 6) G R-^ represents a two-dimensional orbit of stationary solutions u(t) = e^*^^*'^*0(Lij) 
(their time evolution will be specified later) and z{t) G C"'^(M+,Z^) solves the time-evolution 
equation in the form 

ii = {H- uj)z -{9- uj){ct){uj) + z) - iud^(t){uj) + N(0(a;) + z) - N(0(w)), (5) 

where H = —A + V, [N('0)]„ = 7|'(/'„pP^/;n, and duj4>{'^) exists thanks to Lemma [H The 
linearized time evolution at the stationary solution 4>{ijo) involves operators 

L^ = H -uj + W, L+ = H-uj + {2p + 1)W, 

where Wn = ^(thi{tjj) and W decays exponentially as \n\ — > oo thanks to Lemma [TJ The 
linearized time evolution in variables v = Re(z) and w = Im(z) involves a symplectic structure 
which can be characterized by the non-self-adjoint eigenvalue problem 

L+v = —Aw, L_w = Av. (6) 

Using Lemma [H we derive the following result. 

Lemma 3 (Double null subspace). For any e G (0,eo), the linearized eigenvalue problem ^ 
admits a double zero eigenvalue with a one- dimensional kernel, isolated from the rest of the 
spectrum. The generalized kernel is spanned by vectors {0,cj){uj)), {—du,cj){io),0) G l^ satisfying 

L^4>{uj) = 0, L+d^^4>{ijj) = 4>{lo). 
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If (v, w) G is symplectically orthogonal to the double subspace of the generalized kernel, then 



(v,0(a;)) = O, (w,9^0(a;)) = 0, 

where (u,v) := ^nez'^nWn- 

Proof. By Lemma 1 in [T8], operator H has the essential spectrum on [0,4]. Because of the 
exponential decay of W as |n| oo, the essential spectrum of L+ and -L_ is shifted by 
—CO « —ujo > 0, so that the zero point in the spectrum of the linearized eigenvalue problem ([6]) 
is isolated from the continuous spectrum and other isolated eigenvalues. The geometric kernel 
of the linearized operator L = diag(L-|-, L_) is one-dimensional for e € (0, eg) since L-(f){u;) = 
is nothing but the stationary DNLS equation ([2]) whereas L+ has an empty kernel thanks to 
the perturbation theory and Lemma [TJ Indeed, for a small e G (0, eo), we have 

By the perturbation theory, a simple zero eigenvalue of for e = becomes a positive 
eigenvalue for e > (if 7 = +1). The second (generalized) eigenvector {—duicj>{u),0) is found 
by direct computation thanks to Lemma [TJ It remains to show that the third (generalized) 
eigenvector does not exist. If it does, it would satisfy the equation 

L_wo = -d^4>{u:). 

However, 

1 d 

(0(a;),a.(/,(a;)) = - — ||</.(^)||f. = ^^-^ (1 + 0{e)) + 

for e S (0, eo) by Lemma [TJ Therefore, no wq G exists. □ 

To determine the time evolution of varying parameters (w, 0) in the evolution equation ([2|), 
we shall add the condition that z(f) is symplectically orthogonal to the two-dimensional null 
subspace of the linearized problem ([6|). To normalize the eigenvectors uniquely, we set 

||0(u;)|b.' ^2 

and require that 

(Rez(t),^i) = (Imz(t),^2) = 0. (8) 

By Lemma [TJ both eigenvectors i/?]^ and 1/^2 ^.re locally close to "j/jq! the eigenvector of K for 
eigenvalue wq, in any norm, e.g. 

11^1 -^oll/2 + 11^2 -V'oIIp < (^e, (9) 

for some C > 0. Although the vector field of the time evolution problem ([5]) does not lie in 
the orthogonal complement of ■j/jq) that is in the absolutely continuous spectrum of i?, the 
difference is small for small e > 0. We shall prove that the conditions ([8]) define a unique 
decomposition 
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Lemma 4 (Decomposition). Fix e > and 5 > be sufficiently small. Assume that there 
exists T = T[e,S) and Co > 0, such that u{t) G C^([0, T],/^) satisfies 



|u(t) -0(a;o + e))||/2 < Co6e^p , 



(10) 



uniformly on [0,T]. There exists a unique choice of {uj, 6) G C^{[0,T],M?) andz{t) £ C^{[0,T],P) 
in the decomposition ^ provided the constraints (S^ are met. Moreover, there exists C > 
such that 

\uj{t) -ujo-e\ <C6e, \9{t)\<C6, \\z{t)\\i2 < C5e^ , (11) 
uniformly on [0,T]. 



Proof. We write the decomposition dH) in the form 

z = e'^ (u - cl,{uQ + e)) + (e'Ui^o + e) - (t>{^] 



(12) 



First, we show that the constraints ([8|) give unique values of (lo, 9) satisfying bounds (jlip 
uniformly in [0, T] provided the bound (jlOp holds. To do so, we rewrite ([8]) and (jl2p as a 
fixed-point problem F(iv, 6) = 0, where F : i— > is given by 



F(a;,( 



(Re(u-0(°))e*^Vi) + (0^°^ cose-0(w),-0i) 
(Im(u - 0(°))e^^ ^2) + (<^^°^ sin 9, -02) 



where 0^'''' := 4'{ijjQ+e). We note that F is in (0, a;) thanks to Lemma[TJ Direct computations 
give the vector field 

' (Re(u-0(o)),0f^\ ' 



F(u;o + e,0) 

and the Jacobian D'F{u)q + e, 0) = Di + D2 with 
Di = 



(Im(u- ^W),^ 



{0)\ 



Do 



(Re(u - 0(0)), 5^05°)) -(Im(u - 0(°)),0f 

(0 
2 



(Im(u - 0(°)),5^0f ) (Re(u - 0(°)),'0f^^ 



where -j/j^'^g = 0i2U=wo+<: ^'^d duji/^^i2 = 9^^-01 2 U=wo+<:- Thanks to the bound ([T0|) and the 
normalization of 0^ 2' there exists an (e, (5)-independent constant Cq > such that 

||F((Jo + e,0)|| < CoSe^. 
On the other hand, DF{ujo + e, 0) is invertible for small e > since 



.(0) 



|(Di)ii| > Cie^-\ |(Di)22| > C2e^ 



and 



|(D2)ii| + |(D2)2i| < C35e^ \ |(D2)i2| + |(D2)22| < C^Se^, 
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for some (e,(5)-independent constants Ci,C2,C3,C4 > 0. By the Implicit Function Theorem, 
there exists a unique root of F(cj, 6) = near {ljq + e, 0) for any u(t) satisfying (fTO]l such that 

\Lo{t) - iOQ - e\ < C6e, \e{t)\ < C6, 

for some C > 0. Moreover, if u(t) e C^{[0,T],P), then {uj,e) G Ci([0, T], M^). Finally, 
existence of a unique z(t) and the bound ||z(t)||;2 < C6e^p follow from the representation ([T2D 
and the triangle inequality. □ 

Assuming {to, 6) e C^{[0,T],M.'^) at least locally in time and using Lemma [H we define 
the time evolution of {to, 9) from the projections of the time evolution equation ^ with the 
symplectic orthogonality conditions ([8]). The resulting system is written in the matrix-vector 
form 



A(oj, z) 



(13) 



where 



A(lj, z) 



{d^ct){uj),ip-^) - {Rez,d^tpi) (Imz,Vi) 
(Imz, 5^V2) iH'^) + Rez, 



and 



f{uj,z) 



(ImN(0 + z) - Wz,^/7i) 
(ReN(V + z) - N(0) - {2p + l)Wz, ^2) 

Using an elementary property for power functions 



\\a + 6|2P(a + 6) - \a\'^Pa\ < Cp{\a\'^P\b\ + 161^^'+^), 
for some Cp > 0, where a, 6 S C are arbitrary, we bound the vector fields of ([5]) and (I13p by 

||N(0(a;) + z)-N(0H||p < C (|||<^H|2P|z| + ||z||fr^) , (14) 

2 

||f(^,z)|| < Cj;(|||<^(a;)|2p-i|^^.||z|2||,. + |||^^.||zp^'+i||,,), (15) 
for some C > 0, where the pointwise multiplication of vectors on Z is understood in the sense 

(|<^IIV'l)n = <l)n1pn- 

By Lemmas [J and HI A(a;, z) is invertible for a small z G and a small e G (0, eo) so that 
solutions of system satisfy the estimates 



a; 



< Ce'--r{M,\\z\'h + m2\\^\'y) 



|^_^| < Ce'-p{M,\\z\^y + \\\iP^\\z\^\\ii) 
1 

for some C > uniformly in ||z||p < Coe'^p for some Cq > 0. 



(16) 
(17) 
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J_ 

Remark 4. The estimates (1161) and (1171) show that if ||z||;2 < C5e^p for some C > 0, then 



|a;(t) -u;(0)| < CS'e', 



t 

2,2 



9{t) - / L0{t')dt' 







uniformly on [0, T] for any fixed T > 0. These bounds are smaller than bounds (jlip of Lemma 

Co 
<5e 



m They become comparable with bounds (jlip for larger time intervals [0,T], where T < — 



for some Cq > 0. Our main task is to extend these bounds globally to T = oo. 

By the theorem on orbital stability in [24], the trajectory of the DNLS equation ([T|) orig- 
inating from a point in a local neighborhood of the stationary solution (f){io{0)) remains in a 
local neighborhood of the stationary solution 4>{u){t)) for all t G M+. By a definition of orbital 
stability, for any /^q > there exists a i^o > such that if |a;(0) — wqI < then \uj{t) < /Uq 
uniformly on t G R_|_. Therefore, there exists a (5(e) for each e G (0, eo) such that T(e, 5) = oo for 
any (5 G (0, (5(e)) in LemmaSl To prove the main result on asymptotic stability, we need to show 
that the trajectory approaches to the stationary solution cf){uJoo) for some uj oo G (wo,wo + eo)- 
Our main result is formulated as follows. 

Theorem 1 (Asymptotic stability in the energy space). Assume (V1)-(V3), fix j = +1 and 
p > 3. Fix e > and 6 > be sufficiently small and assume that 6{0) = 0, ti;(0) = loq + e, and 

||u(0) -0(u;o + e)||/2 < CoSe^ 

for some Cq > 0. Then, there exist Wqo G {ujo,ujo + eo), {uj,6) G C^(M+,M^), and a solution 
u{t) eX := Ci(M+,/2) n L^(M+,Z°°) to the DNLS equation (C]) such that 

lim Lo{t) = Woo, ||u(i) - e-^^W0(a;(t))||x < C(5e^/(2p)_ 

Theorem [T] is proved in Section 4. To bound solutions of the time-evolution problem ([SJ 
in the space X (intersected with some other spaces of technical nature), we need some linear 
estimates, which are described in Section 3. 

3 Linear estimates 

We need several types of linear estimates, each is designed to control different nonlinear terms 
of the vector field of the evolution equation ([5]). For notational convenience, we shall use 
and In to denote space on t G [0, T] and I'' space on n G Z, where T > is an arbitrary time 
including T = oo. The notation < n >= (1 -|- n^)^/^ is used for the weights in In norms. The 
constant C > is a generic constant, which may change from one line to another line. 

3.1 Decay and Strichartz estimates 

Under assumptions (V1)-(V2) on the potential, the following result was proved in [18j. 

Lemma 5 (Dispersive decay estimates). Fix a > | and assume (V1)-(V2). There exists a 
constant C > depending on V such that 

\\{n)-'^e-^'"PaAH)A\i2 < C(l + t)-3/2||(^)-f|| (18) 

II 1 1 1,,^ n 

||e-**^P,.e.(^)f|Loo < C(l + t)-V3||f||i (19) 
for all t G M+, where Pa.c.{H) is the projection to the absolutely continuous spectrum of H. 



8 



Remark 5. Unlike the continuous case, the upper bound U9\) is non-singular as t ^ because 
the discrete case always enjoys an estimate ||f||;oo < ll^lb^ < ll^lb^- 

Using Lemma[5]and Theorem 1.2 of Keel-Tao [lOJ, the following corollary transfers pointwise 
decay estimates into Strichartz estimates. 

Corollary 1 (Discrete Strichartz estimates). There exists a constant C > such that 

'|e-^*^P„.,.(F)f||^e,oonLoop < Cm\i2, (20) 



f e-"^'-'^'' Pa.c.{H)g{s)ds 
Jo 



where the norm in L^ln is defined by 



< C\\ghu2, 



1/p 



(21) 



3.2 Time averaged estimates 

To control the evolution of the varying parameters {uj,9), we derive additional time averaged 
estimates. Similar to the continuous case, these estimates are only needed in one dimension, 
because the time decay provided by the Strichartz estimates is insufficient to guarantee time 
integrability ofio{t) and 9{t)—uj{t) bounded from above by the estimates (fTBll and (fTTll . Without 
the time integrability of these quantities, the arguments on the decay of various norms of z(t) 
satisfying the time evolution problem ([5]) cannot be closed. 

Lemma 6. Fix o" > | and assume (VI) and (V2). There exists a constant C > depending 
on V such that 



< n >-3/2 e-^*^P„.a/7)f ll;^,^. < C\ 



itH 



Pa.c.iH)Fis)dt 



<n> 



<n> 



I e-'^'-'^'' PaAH)F{s)ds 
Jo 

f\-^(t-s)Hp^^^H)¥{s)ds 
Jo 

f e-"^'~'^''Pa.cXH)Fis)ds 
Jo 



loo J 2 



Lfl?^nL^ll 



< c\ 

< c\ 

< c\ 

< c\ 



< n >2/2 F||,w2, 



< n F||,i 

< n >^ F 



;1 r 2 



(22) 
(23) 

(24) 

(25) 

(26) 



To proceed with the proof, let us set up a few notations. First, introduce the perturbed 
resolvent Rv{\) '■= {H — X)~^ for A G C\[0, 4]. We proved in [18| Theorem 1] that for any fixed 
oj G (0, 4), there exists Ry{oj) = linieio R{'^ =t ^e) in the norm of B{a^ —a) for any o" > ^, where 
B{a, —a) denotes the space of bounded operators from to 

Next, we recall the Cauchy formula for e**^ 

e-'^^ Pa.c.{H) = - [ e-'^lmRv{uj)duj = — [ e"**"^ [i?+(w) - R~{uj)] du, (27) 
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where the integral is understood in norm B(a,—a). We shall parameterize the interval [0,4] 
by a; = 2 - 2cos(6') for e [-7r,7r]. 

Let XOj X S Co° : Xo + X = 1 for all 9 G [— tt, vr], so that 



suppxo C [-6*0, 6*0] U (-vr, -TT + 6*0) U (vr - ^o, ^r) 



and 



suppx C [60/2, vr - 60/2] U[-7T + 9o/2, -9o/2] , 

where < ^0 ^ f • Note that the support of x stays away from both and vr. Following 
Mizumachi [T5] , the proof of Lemma E] relies on the technical lemma. 



(28) 
(29) 



Lemma 7. Assume (VI) and (V2). There exists a constant C > such that 

SUp||x«y(w)f||L2(o,4) <C||f||,2, 

sup II < n >-3/2 x^R^ii^MLliOA) < Cmiii. 

The proof of Lemma [7| is developed in Appendix A. Using Lemma U\ we can now prove 
Lemma [6l 

Proof of Lemma\B Let us first show (p^ . since it can be deduced from p^ . although, it can 
also be viewed (and proved) as a dual of (|22|) as well. Indeed, (|24p is equivalent to 



<n>' 



J 



By the Krein's theorem, for every Banach space X, the elements of the space /^(X) are weak 
limits of linear combinations of functions in the form (5nno^) where x € X, ng G Z, and 



is Kronecker's symbol. Thus, to prove the last estimate, we need to check if it holds for 



Gnis) = Sn,no9is), where g £ Lf. By Minkowski's inequality, the obvious embedding / 
and the dispersive decay estimate (fTSl) for any a > |, we have 



<n> 



< C 

< C 



Jo 

Jo 

<C\\g\\L., 



/2 



\g{s)\ds 



\gis)\ds 



{l + t- s)3/2 



where in the last step, we have used Hausdorff- Young's inequality * ^ L^. 

We show next that ([23]), ([25]), ([26]) follow from ([22]). Indeed, ([23]) is simply a dual of ([22]) 
and ([23]) is hence equivalent to ([22|) . For ([25]) . we apply the so-called averaging principle, which 
tells us that to prove ()25p . it is sufficient to show it for F{t) = 6{t — to)¥, where f G and 
^{i ~ to) is Dirac's delta-function. Therefore, we obtain 
t 



<n> 



-i(t-s)H 



6{s-to)Pa.cXH)ids 



<n> 



loo T 2 

'"n 



e-^(*-*«)^P,.,,(i^)f|| 



< II < n >-=^/2 e-'^^^-^o^^ Pa,c.{H)i\\ 
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where in the last step, we have used ()22| 
For ()26p . we argue as follows. Define 



TF{t) = [ e-'^'-'^'' Pa.c.{H)F{s)ds 



-itH 



-isH 



where f = f^e~^^^ Pa.c.{H)^{^)'^^- By an application of the Strichartz estimate ()20l) and 



subsequently (j23j) . we obtain 

WTFWL^ii^nL^il < C||f|b^ < II < ^ n^L^ <C\\<n >^ F||,.^2 = C\\ < n >^ F||^2,. , 

where in the last two steps, we have used Holder's inequality and the fact that and 
commute. Now, by the Christ-Kiselev's lemma (e.g. Theorem 1.2 in [IQj), we conclude that 
the estimate ([MD applies to /g e-*(*-'')^Pa.c.(i^)F(s)(is, similar to TF{t). To complete the 
proof of Lemma [71 it only remains to prove (|22p . Let us write 

e-'''' Pa.cXH) = xe-^'^'PaAH) + Xoe-^*^Pa.c.(^) 
Take a test function g(t) such that ||g||;i/,2 = 1 and obtain 

\{xe-''''Pa.c.iH)i,git))n,t\ = - f\xl^Rv{^)i, [ e-''^git)dt)ndu; 

< C [ {\xRv(.uj)i\,\g{u^)\)ndu; 

Jo 

< 1 1 X^y ('^) f 1 1 «5f L2 (0,4) 1 1 g I bi L2 (0,4) • 

By Plancherel's theorem, ||g||/iL2 (0,4) < l|g|biL2(R) < llslLiL? = 1. Using ([28]), we obtain 
||xe-**^Pa.c.(^)f|Lc»i2 = sup |(xe-^*^Pa.c.(^)f,g(i))n,t| < Cm\i2. 

* l|g|l;l 1,2=1 



Similarly, using (i29]l instead of (128]) . one concludes 

sup |(xoe-**^Pa.c.Wf,g(t))n,t| < C||f|b2. 



<n>-3/2 xoe-''''Pa.c.{H)i 



<n>3/2g|| 1 2=1 



Combining the two estimates, we obtain (|22p . 

4 Proof of Theorem [1] 

Let y(t) = e~''^^^^z{t) and write the time-evolution problem for y{t) in the form 

iy = ffy + gi + g2 + g3, 
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where 

gi = (n(<^ + ye-'^) - N(0)) e-^^ g2 = -(^ - u;)<^e-'^ gg = -id;dM^)e-'' ■ 



Let Pq = (•,'0o)'0oi Q = ~ Pq) = Pa.cX^)^ s-^d decompose the solution y(t) into two 
orthogonal parts 

y(t) =a(t)^o + ^(0, 

where {tpQ,r]) = and a{t) = (y(t), ■j/jq). The new coordinates a{t) and rj{t) satisfy the time 
evolution problem 

id = woa + (g,V'o)' 
irj = Hrj + Qg 

where g = Yl^=iSj- The time-evolution problem for rj = Pa.c.{H)ri can be rewritten in the 
integral form as 

rj{t) = e-**^QT7(0) - i [ e-^(*-^)^Qg(s)ds, (30) 
Fix o" > I and introduce the norms 



Mi = \\v\\lIi^, ^2 = NIIl-/2, M3 = \\ <n> riWi^Lh 



^4 = ||a||r2, M5 = ||a||L?o, Mq = \\UJ - UJ{U)\\L 



where the integration in is performed on an interval [0, T] for any T E (0, oo). Our goal is 
to show that to and 6 — oj are in L\, while the norms above satisfy an estimate of the form 



Y,M,<C\\ymil+c[Y,MA (31) 

and 

Me < Ce^~p{M3 + M^f , (32) 

for some T-independent constant C > uniformly in ^ Mj < CSe^p , where small positive 

i=i 

j_ 

values of (e, 5) are fixed by the initial conditions w(0) = ujq + e and ||y(0)||i2 < Co^e^p for 
some Co > 0. The estimate ([3T]) and (j32|) allow us to conclude, by elementary continuation 
arguments, that 

5 

Y,^I,<C\\Ym\il<C5e^- 
i=i 

and \ijj{t) — ujQ — €\ < Cb'^e^ uniformly on [0, T] for any T G (0, oo). By interpolation, a G so 
that z(t) G -^^^([0, T], V^). Theorem [1] then holds for T = oo. In particular, since G L|(M+) 
and |a;(i) — loq — e\ < C5^e^, there exists Woo := limf_>oo <^(i) so that Woo £ ('^Oi'^o + co)- In 
addition, since z{t) G L^(M+,^5f ), then 

hm ||u(t) - e-^^W0(a;(t))|bc. = lim ||z(t)||;c>o = 0. 
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Estimates for Mq: By the estimate (fT6]) . we have 

/ \dj\dt < Ce^'p\\ Kny-^" \y\'^\\i^Li {\\ Kny^" -tp^Wii + \\ <n>^^ tP^h) 
Jo " ' 

< Ce^-p\\<n>-'' y\\%^2 

'■n 

< Ce^"p(M3 + M4)2, 

where we have used the fact that ij^i and i/'2 decay exponentially as |n| oo. As a result, we 
obtain 

Me < 11^11^1 < Ce^~p{M3 + M^f. 

Similarly, we also obtain that 

/ \6-uj\dt<Ce^'p{Ms + Mif. 
Jo 

Estimates for M4 and M5: We use the projection formula a = (y, i/;q) and recall the 
orthogonality relation ([8]), so that 

(z, -00) = (Rez, -00 - -01) + i(Imz, Vo - ^2)- 

By Lemma [1] and definitions of i/'i 2 i'^ ©> have 

II < n >2'^ (^0 - V'i,2)lb2 < C\u; - uJo\ 

for some C > 0. Provided cr > |, we obtain 

^4 = |Ky,^o)llL? < IKRez,V'o-^i)llL? + liamz,V'o-V^2)llL? 

< II < n >-2" z||i2;2 (11 < n (Vo - V'i)llL-i2 + II < n (t/^g - ^2)llLr'0 

< C\\ < n y-'' y\\ic^^2\\uj - loqWl^ <C{M'i + Ma)Mq 
and, similarly, 

< llyllir'Si (iK^o - ^i)llL-i2 + 11(^0 - ^2)\\Lri^}) < ciAh + M5)Mq. 

Estimates for M3: The free solution in the integral equation (j3U|) is estimated by (|22|) as 
II < n >-'^ e-^*^Qr7(0) 11,5^3^2 < || < n >-3/2 6^*^^77(0)11,5.^2 < C||T,(0)|b2. 

Since lj and 6 — u are thanks to the estimates above, we treat the terms of the integral 
equation ([30]) with g2 and gs similarly. By (j25|) . we obtain 

II <'^>-'' re"*^*-'^''Qg2,3(sM5||/5J=L2 < C||g2,3||ij;2 
t/ 

< C (^11^ -a;||^i||^(w)||ic«,2 + 11^11^1 ||a^0(u;)||ic«,2) 

< Ce^"^(M3 + M4)2. 

13 



On the other hand, using the bound (|14p on the vector field gi , we estimate by (j24p and (p5 



'i{t-s)H 



<c(\\<n yWi^^.W <n>^ |0(u;)P^||l^,i + llallj+i llV-oll'li+i) + Nll'Ci,^)) 
< C ({Ms + M^)M, + MiM?/-' + WvfXl., 



2{2p+l) I ) 



where we have 



and 



|2p+l 



|2P-1||„||2 



2p+l ^ II'^IIloo Il'^llr2- 



<n>'' |(/)(tj)|2P||,i < C\\uj -LOo\ 



the latter estimate follows from Lemma [TJ 

To deal with the last term in the estimate, we use the Gagliardo-Nirenberg inequality, that 
is, for all 2 < r, li; < cxD such that - + — < 1, there is a C > such that 



WvUnr, < c [\\v\\L<ii^ + Wvhri?, ) = c{Mi + M2). 

If p > 3, then {{2p + 1), 2(2p + 1)) is a Strichartz pair satisfying 2^fi + 2{2p+i) — ^ ^^'-^ hence, 
combining all previous inequalities, we have 

M3 < C {\\v{0)\\i2 + e^^p (M3 + M^f + (M3 + M4)M6 + mIm^P'^ + (Mi + M2)=^p+^) , 

which agrees with the estimate (|3ip for any p > 3. 

Estimates for Mi and M2: With the help of (j20p . the free solution is estimated by 

\\e-''''Qv{0)\\Lti?^nLril<CMO)\k. 
With the help of (j2ip . the nonlinear terms involving g2,3 are estimated by 



Qg2,3is)ds 



< C'l|g2,3|lLi«2 



< Ce^"i(M3 + M4)2. 



The nonlinear term involving gi is estimated by the sum of two computations thanks to the 
bound ()14p . The first computation is completed with the help of (I26p . 



-i{t-s)H 



Q|0(u;)|2P|y|ds 



Lfl?pnL^ll 



< C||<n>3|c/,H|2P|y||| , 



< II < n >3+- |(^H|^^||ic»,2 II < n y||,^i2 

< CiMs + M4)Me, 



|2P| 
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whereas the second computation is completed with the help of ()2ip 



^_,(t_,)Hg|y|2p+lj^ 



< 



C7|||y|2p+i 




< c [mImI^- 



+ (Afi + Af2)'^+^), 



provided p > 3 holds. We conclude that the estimates for Mi and M2 are the same as the one 
for M3. 

5 Numerical results 

We now add some numerical computations which illustrate the asymptotic stability result 
of Theorem [TJ In particular, we shall obtain numerically the rate, at which the localized 
perturbations approach to the asymptotic state of the small discrete soliton. One advantage of 
numerical computations is that they are not limited to the case of p > 3 (which is the realm 
of our theoretical analysis above), but can be extended to arbitrary p > 1. In what follows, we 
illustrate the results for p = 1 (the cubic DNLS), p = 2 (the quintic DNLS), and p = 3 (the 
septic DNLS). 

Let us consider the single-node external potential with Vn = —Sn,o for any n ^ 'L. This 
potential is known (see Appendix A in [12]) to have only one negative eigenvalue at < 0, 
the continuous spectrum at [0,4], and no resonances at and 4, so it satisfies assumptions 
(V1)-(V3). Explicit computations show that the eigenvalue exists at = 2 — \/5 with the 
corresponding eigenvector -00,™ = e"'''"'! for any n G Z, where k = arcsinh(2~^). The stationary 
solutions of the nonlinear difference equation (l2|) exist in a local neighborhood of the ground 
state of H = —A + V, according to Lemma [TJ We shall consider numerically the case 7 = 
— 1, for which the stationary solution bifurcates to the domain iv < loq. Figure [T] illustrates 
the stationary solutions for p = 1 and two different values of to, showcasing its increased 
localization (decreasing width and increasing amplitude), as oj deviates from ujq towards the 
negative domain. 

In order to examine the dynamics of the DNLS equation ([T]) we consider single-node ini- 
tial data Un = A6n,o for any n E Z, with A = 0.75, and observe the temporal dynamics of 
the solution u(t). The resulting dynamics involves the asymptotic relaxation of the localized 
perturbation into a discrete soliton after shedding of some "radiation". This dynamics was 
found to be typical for all values of p = 1,2,3. In Figure [21 upon suitable subtraction of the 
phase dynamics, we illustrate the approach of the wave profile to its asymptotic form in the 
l°° norm. The asymptotic form is obtained by running the numerical simulation for sufficiently 
long times, so that the profile has relaxed to the stationary state. Using a fixed-point algo- 
rithm, we identify the stationary state with the same l'^ norm (as the central portion of the 
lattice) and confirm that the result of further temporal dynamics is essentially identical to the 
stationary state. Subsequently the displayed norm of the deviation from the asymptotic 
profile is computed, appropriately eliminating the phase by using the gauge invariance of the 
DNLS equation i^. 

We have found from Figure [2| in the cases p = 3 (top panel) , p = 2 (middle panel) and 
p = 1 (bottom panel) that the approach to the stationary state follows a power law which is 
well approximated as oc t~^/'^ . The dashed line on all three figures represents such a decay in 
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Figure 1: Two profiles of the stationary solution of the nonlinear difference equation ([2|) for 
Vn = —Sn,o, P = 1, and for uj = —2 (solid line with circles) and lo = —5 (dashed line with stars). 

each of the cases. We note that the decay rate observed in numerical simulations of the DNLS 
equation ([1]) is faster than the decay rate oc for any p > in Theorem [TJ 

A Proof of Lemma [7] 

For the proof of Lemma [71 we will have to show both the "high frequency" estimate (j'28|) and 
the "low frequency" estimate (|29p . To simplify notations, we drop the bold- face font for vectors 
on Z in the appendix. 

A.l Proof of (EH]) 

Recall the finite Born series representation of Ry 

Riu) = Roiu) - Ro{u;)VRoiu;) + Ro{u;)V Riio)V Roico) , (33) 

which is basically nothing but the resolvent identity iterated twice. We have shown in |18j that 
for the "sandwiched resolvent" Gu,w{^) = URv{u^)W , we have the bounds (see estimate (33) 
in [Ml) 

<C\\U\\i2\W\\i2. (34) 

for any fixed cr > |, where u> = 2 — 2cos(0). 

For the three pieces arising from (|33p . similar arguments apply. Starting with the free 



sup 

6»ef-7r,7 



dO 
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Figure 2: Evolution for p = 3 (top), 2 (middle), 1 (bottom) of ||u(t)— e~*^(*^(/)(u;oo)|| as a function 
of time in a log-log scale (solid) and comparison with a t^'^/^ power law decay (dashed) as a 
guide to the eye. 



resolvent term, we have 



< C sup 



Introducing the sequence 



sup / xliRoi^^) f)n\'^duj < Csup / 

I 



X 



sin(0) 



^id\m-n\ j!^ 



mG2 



d0 < 







2 






( 




+ 




dO 




m>n 




m<ra 





M ■- 



fm rn>n 
m < n 

|2 



we see that the last expression is simply C( lb" 1 1^2 [51^/2 n-eo/2] 
is equal by Plancherel's identity to 



+ 



Cml + \\f-gm<2C\\f\\l. 



For the second piece in (133]) . we use that ||ii^(tj)||^i^;oo < C/ sin(^) and |sin(^)| > Co on 
[9o/2,TT — 9o/2] for some Co > 0, to conlcude 



sup\\xR^iu;)VR^iu;)f\\l,,,^ 



< 



X 



n sin^(^) 



nfc£, 7 — TT 



< C||y|biSUp/ x|(^?M/)nrd^, 



by the triangle inequality. At this point, we have reduced the estimate to the previous case, 
provided that V ^ l^. 

For the third piece in we make use of We have, similar to the previous estimate, 
sup \\xR^{o,)VR${Lo)VR^{u)f\\l,^,^,) = 
sup \\xR^iu;)VR^{u;)\V\'/hgn{V)V'/^R^{u;)f\\l,^,^,^ = 
sup||xiio^(u;)G,.,H,|V2,g,(v')[|y|i/2i?o^(u;)/]||i.(o_4) 

<C^ll^ll|lll^l'/'ll|lll^l'/'ll'iSup r x\{R^iu;)f)n\^de, 

where in the last inequality, we have again reduced the estimate to the first case. 

A.2 Proof of am 



We only consider the interval [—9o,Oq] in the compact support of Xo(^) since the arguments 
for other intervals are similar. Following the algorithm in [15j and the formalism in [18], we let 
'4>^{0) be two linearly independent solutions of 



V'n+l + V'n-l + (W - 2)'4jri = KV'n, n G Z, 



(35) 
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according to the boundary conditions \ipn — e^*"^| ^ as n ^ ±00. Let V'n (^) = e^*"^^^(0) 
for all n € Z. Using the Green function representation, we obtain 



n 



The discrete Green function for the resolvent operators R^{uj) has the kernel 



[R 



V 'J n,m 



V'^(6'±)V'm(^±) n>m 

W{d±) \ iJ+iO±)^niO±) for n < m 



where 9- = -9+, 9^ G [0,7r] for uj G [0,4], and W{9) = V"] = V-^C+i - ^n+iV'n is the 

discrete Wronskian, which is independent of n G Z. We need to estimate 



\\xoRH^)f\\l2^oA) 



^ 2x1 sin 9d9 



n-l 



J2 i^imMU+Yl ^n{e)&)U 



We may assume that n > 1 for definiteness and split 

ra-l n-l -1 -1 

i^-{9)fm=Y.i';niO)fm+ Yl ^''^'fra+ e'^' - 1) fm := h + h + h 



m=—oo 



m=0 



m=—oo 



m=—oo 



and 



We are using the scattering theory from [TB] to claim that 

sup iW^^mii^iz^) + \\{n)-'^^m\i^iz^)) < oo, 
ee[-eo,eo] 



(36) 



where (n) = (1 + "n?)^^"^. Then, we have 



'n-l 



1/2 



1 



1/2 



\m=0 / 

i/3i < f E 



EI/-I 

\rn=0 
\ 1/2 



<Ci(n)3/2||/||p, 



E I/" 



1/2 



>,m=— oo 



141 < Ei^-(^)-ii 



1/2 



>,m=— oo 



El/" 



1/2 



< a 



k=—oo 

oo 

Y\k-m\\Vk\ 



l=m 



for some Ci, C3, C5 > 0. We note that 



E \^-k\\^k\ 



k=—oo 



< 



A;=— 00 



< G 



411 ^ ll/i, 
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for some C4 > 0. Therefore, the brackets in and Is are bounded if y € l^^ for a > |. Since 
I2 and I4 are given by the discrete Fourier transform, Parseval's equahty imphes that 

J —TT 

for some C2 > 0. Using now the fact that > Wq and | sin9\ < Cq uniformly in [—9o, Oq], 

the support of Xo(^)i and using the property ([36]) . we obtain 

||xoi?^M/||i. (0,4) <C{1 + {nf + {nf) 

which gives ()29p . 
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